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Abstract. We present various aspects of the saturation model which provides good 
description of inclusive and diffractive DIS at small x. The model uses parton 
saturation ideas to take into account unitarity requirements. A new scaling predicted 
by the model in the small x domain is successfully confronted with the data. 



1. Introduction 

Deeply inelastic electron-proton scattering (DIS) at small value of the Bjorken variable x 
(< 1) attracted a lot of attention, mostly due to the experimental results from HERA. 
From the theory side, the small-x DIS opens a new kinematic regime for the QCD 
studies. The predicted by QCD, and confirmed by the data |IJ, strong rise of the 
proton structure functions with decreasing x is the indication that at small x the proton 
structure is dominated by dense gluon systems. In DIS at moderate values of x, the linear 
QCD evolution equations lead to good description of this process, explaining scaling 
violation. At small x, however, the problem is more complicated since recombination 
processes between partons in a dense system have to be taken into account. At the 
formal level, these processes allow to restore unitarity of the DIS cross sections as x — > 0, 
violated in the QCD description based on the linear evolution equations, by taming the 
rise of the parton distributions. This effect is called parton saturation [[| . 

It was unclear for a long time whether recombination processes (or unitarization 
procedures in general) are important in the HERA kinematic range. Most of the 
analyses of this problem were based on the results obtained in the double logarithmic 
approximation in which x is small but Q 2 is large. The linear evolution equations in 
this case are modified by non-linear terms describing parton recombination. In such 
approximation it was found that the impact of non-linearity is either small or masked 
by the choice of initial conditions for the linear evolution j3|. 

The breakthrough came from realization that the proper starting point is the 
high energy QCD factorization formula for the structure functions, derived under the 
assumption that x is small but Q 2 is arbitrary (provided Q 2 ^> A 2 and pQCD is 
justified). This allows to find the necessity for unitarity, realized in the spirit of parton 
saturation, in the transition from DIS to small Q 2 scattering M. The full potential of 
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the unitary description is revealed for diffractive DIS ||. In particular, the constant 
ratio of diffractive to total DIS cross sections as a function of x and Q 2 is naturally 
explained. 

The approach we are going to present is partially phenomenological. We take into 
account several important features of parton saturation in the parameterization of the 
QCD interactions in a dense gluon system, especially the existence of a saturation scale. 
The proposed parameterization is very efficient in the description of the data and may 
serve as a guidance for the detailed QCD studies. 

2. Saturation model 

The saturation model was formulated and compared at length to DIS data in |4], ||. 
Here we describe this model, discussing some details which were not presented in the 
original formulation. For related approaches see @. 

In the rest frame of the proton, the QCD description of DIS at small x can be 
interpreted as a two-step process. The virtual photon (emitted by the incident electron) 
splits into a qq dipole which subsequently interacts with the proton. The proton 
structure function F<i results from the high energy QCD factorization theorem. In 
terms of virtual photon-proton cross sections <Jt,l f° r the transverse and longitudinal 
polarized photons 

F 2 (x,Q 2 ) = Q 2 /4Tr 2 a em (a T + a L ) (1) 

and 

o~t,l = J d 2 rdz\q T<L (r,z,Q 2 )\ 2 &{x,r), (2) 

where ^t,l is the light-cone wave function of the virtual photon and a is the dipole 
cross section describing the interaction of the qq dipole with the proton. In equation (|2|) 
r is the transverse separation of the qq pair and z is the photon's momentum fraction 
carried by the quark, see Figure [l|. Thus (r, z) are good quantum numbers conserved 
by the interaction in the considered approximation. 




Figure 1. Schematic representation of the basic factorization in inclusive DIS at 
small x. 
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The photon wave functions are known from pQCD, e.g. for transverse photons and 
massless quarks 



\*t{t,z,Q 



2\|2 
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zf\ Q Kf(Qr), 



(3) 



where Q = z(l — z) Q 2 and Ki is the Bessel function. The computation of the 
dipole cross section a has been attempted within pQCD assuming different types of 
net colourless gluon exchange (e.g. DGLAP or BFKL ladders). Most of these attempts, 
however, are plagued by problems with unitarity of finally computed cross sections. 
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Figure 2. Dipole cross section for different values of x. 



In our approach we built in unitarity in the dipole cross section by proposing the 
following phenomenological form 

a(x,r) = ao{l-exp(- i ^- y )}, (4) 

with Ro(x), called saturation radius, given by 

where Qq = 1 GeV 2 . The parameters o"o = 23 mb,xo = 3 • 10~ 4 and A = 0.29 were 
found from the fit to all inclusive DIS data at x < 0.01. At small r, a features colour 
transparency, a ~ r 2 , which is purely pQCD phenomenon. For large r, saturation 
occurs, a ~ o~q. The fact that a is limited by the energy independent cross section may 
be regarded as a unitarity bound. The transition between the two regimes is governed 
by Ro{x). As illustrated in Figure 0, for x — > the transition occurs for decreasing 
transverse sizes. This is an essential feature of parton saturation pioneered by the 
analysis [f2j in which internal saturation scale Q s {x) appears. We built in such a scale 
into our model: Q s (x) ~ l/Ro(x). 

With the proposed dipole cross section we achieved good description of the DIS 
data at small x, including the transition to small Q 2 values, see Figure |3|. The reason 
for this can be easily understood after performing the following qualitative analysis. 
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Figure 3. ct 7 * p = <tt + <Jl for different energies W of the 7*p system. The solid lines 
show the fit results with a light quark mass m/ = 140 MeV while the dotted lines with 
m,f — 0. The line across the curves indicates the position of the critical line. 



3. Qualitative analysis 

We concentrate on the transverse cross section o~t which dominates in F2. The dominant 
contribution to a? results from the behaviour of the Bessel function in (0): K\{x) — 1/x 
for x <C 1 while for x 3> 1, K\{x) is exponentially suppressed. Thus using (|3]), 

a T ~ / — dz[z 2 + (l- z) 2 } &(x, r) 6 [^(1 - z) Q 2 r 2 < l] , (6) 
Jo r Jo 

where Q(x < 1) equals 1 if x < 1, or otherwise. For < r < 2/Q, the theta 
function does not impose any restriction on z and the z-integration gives a constant. 
Such a configuration, called symmetric, is rather uniform in z with the mean value 1/2. 
For large transverse separations, r ^> 2/Q, the theta function heavily restricts z to 
small values: z < l/(Q 2 r 2 ). Now, the ^-integration gives the factor 2/{Q 2 r 2 ). In this 
configuration, called aligned jet, z or (1 — z) ~ 0. Thus, one of the quarks follows the 
photon direction while the other stays with the proton. 
Finally, we obtain 



a r ~ / ^a(x,r) + y — ) a(x,r), (7) 



symmetric aligned jet 
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where we neglected multiplicative numerical factors. For the forthcoming analysis let 
us approximate (|6|) by 

{a Q r 2 /AR 2 (x) for r < 2R (x) 
(8) 
(T for r > 2R (x) , 

which form contains all essential features of the exact formula. The leading Q 2 -behaviour 
of o"t depends on the relation between two scales: the characteristic size of the qq dipole 
l/Q and the mean transverse distance between partons given by Ro{x). 
If the partonic system is dilute l/Q -C Ro(x), and from (0J|) we find 





2/Q<r<2R 

where we indicated the r-integration regions. In terms of the structure function ([]]) we 
obtain scaling with logarithmic violation. Notice that all dipole sizes, including those 
from the non-perturbative region, contribute to scaling. The assumption that R^ ~ x x 
leads to the power-like behaviour, F 2 ~ x~ x , observed in the data. 

When the system of partons becomes dense for the dipole probe and l/Q 3> Ro(x), 
a different behaviour is found 

a T ~ jJo^ + (7 In (^QTj^j + ■ ( 10 ) 

r<2R * v ' r>2/Q 

2R <r<2/Q 

Now, the structure function is in agreement with unitarity: F 2 ~ Q 2 o~q ln(l/x). The 
change of the behaviour of o~t from (H) to ( p^0|) when Q 2 increases while W 2 stays fixed 
(x = Q 2 /W 2 — > 0) is shown in Figure |3|. The same transition is obtained in the Regge 
limit of DIS, Q 2 fixed and W 2 — > oo, since the j*p center-of-mass energy W enters the 
description only through the Bjorken variable x. 

The transition between the scaling and unitary behaviour of the structure function 
depends on x and occurs in the region of the (x, Q 2 )-plane marked by the critical line: 
l/Q = Ro(x). In our interpretation, at the critical line the characteristic size of the 
dipole probe equals the mean transverse separation between partons. The characteristic 
feature of parton saturation is that the transition occurs at increasing values of Q 2 when 
x — > 0. We found that at HERA the transition region is situated for Q 2 1 — 2 GeV 2 , 
which justifies the use of pQCD. 

4. Photoproduction limit 

It is interesting to consider a formal limit Q 2 — ► in the saturation model. The analyzed 
cross sections are divergent in this limit if m/ = 0. However, if a non-zero quark mass 

'5 



is assumed, the limit can be performed and we find for m 2 f 3> Q 2 — > 0, 
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Figure 4. The effective pomeron intercept ap = 1 + A as a function of Q , found from 
the dependence: <7 7 » p ~ (W 2 ) x . The solid line corresponds to the saturation model 
([|) while the dashed line shows the effect when the model is modified to include the 
DGLAP evolution (Section ^J). The data points are from ZEUS. 



OL ~ cr , (12) 
m f 

where we additionally modify the Bjorken variable formula to allow for the 
photoproduction limit 

Q 2 + 4:m 2 f / 

As expected, the longitudinal cross section vanishes when Q 2 = 0. We also see that uif 
plays a crucial role for the value of the transverse cross sections. In our analysis we set 
uif — 140 MeV to obtain good agreement with the HERA photoproduction data. For 
Q 2 ^> m 2 , the light quark mass does not play a significant role. 

From the discussion in Section |3| we know that the energy dependence of <7t changes 
from o~t ~ (W 2 ) x for large Q 2 to o~t ~ ln(iy 2 ) in the photoproduction limit. It appears 
that for each Q 2 we can effectively parameterize the energy dependence through the 
power-like behaviour: ~ {W 2 ) ap ^ 2 ^~ 1 . The found dependence ctp(Q 2 ) is shown 
in Figure f|. Interestingly, ap(Q 2 ) interpolates between the soft and hard pomeron 
intercept values. 
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5. Geometric scaling 



Let us realize that the dipole cross section (^) is a function of the combination r/Ro(x) 
instead of x and r separately, 

a{x,r) = &{r/Ro{x)). (14) 

This has profound consequences for the total cross section <T 7 . p = o~t + &l, if the quark 
mass m/ in the photon wave functions is neglected. Rescaling the integration variable 
r — > r/Ro in (^|), we find that <r 7 * p becomes a function of the dimensionless variable 
r = Q 2 R^(x), being the ratio of the two geometric scales discussed in Section [| 

a rp (x,Q 2 ) = a rp (Q 2 R 2 Q (x)) . (15) 

The non-zero light quark mass, introduced to extrapolate the model down to the 
photoproduction region, does not lead to a significant breaking of a new scaling (which 
we call geometric scaling) in the small x domain H. So does the charm contribution, 
discussed in detail in ||. In Figure reproduced from ||, we illustrate geometric 
scaling for the small-x data in a broad range of Q 2 . 

In its essence, geometric scaling is a manifestation of the presence of the internal 
saturation scale characterizing a dense partonic system, Q s (x) ~ 1/Rq(x). The presence 
of such scale emerges from a pioneering work of f2|, which was subsequently analyzed 
and generalized in |§. In reference flO| the scaling property similar to that postulated 
in the dipole cross section (f|) was found and analyzed in detail in |TT|. 
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Figure 5. Experimental data on cr 7 » p from the region x < 0.01 plotted versus the 
scaling variable r = Q 2 Rq(x). Q 2 values are between 0.045 and 450 GeV 2 . 
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Figure 6. Diffractive qq and qqg contributions. 



6. DGLAP evolution in the saturation model 

In the limit of large Q 2 values, the structure function ([[J) is dominated by the small size 
transverse contribution. In this limit, the high energy formula ([!]) should make a contact 
with the DGLAP formula for F 2 . This allows to find the following relation between the 



dipole cross section at small r and the ordinary gluon distribution JT2 

7T 2 

a(x, r) ps — r 2 a s xg(x, C/r 2 ). (16) 

The gluon distribution g(x,n 2 ) obeys the DGLAP evolution equations. The parameter 
C in fllB]) , as well as the argument of a s , cannot be determined in the considered leading 
logQ 2 approximation. 

In the saturation model a ~ r 2 / R^{x) at small r. Thus the logarithmic 
dependence on r, resulting from the DGLAP evolution of the gluon, is not included 
in the model. In our recent analysis |TJ| we include the DGLAP evolution by proposing 



the following modification of the saturation model 

v x /, / 7T 2 r 2 a s (n 2 )xg(x,fi 2 ) \ \ 

a(x,r) = o-q j 1 - exp I — 1 j , (17) 

where the scale fi 2 = C jr 2 + /i 2 ,, and the parameters C and /i 2 , are determined from the 
best fit to the DIS data. Additionally, two parameters of the initial gluon distribution: 
xg = A g x~ Xg are fitted. The result of the comparison with the data is shown in Figure § 
as the dotted line. As expected, the proper DGLAP limit of a significantly improves 
agreement at large values of Q 2 without affecting the physics of saturation responsible 
for the transition to small Q 2 . The full discussion of the results will be presented in |13l . 



7. Saturation and DIS diffraction 



DIS diffraction, 7*p — > pX, is a good test of the parton saturation ideas incorporated in 
(|j). In the simplest case the diffractive system X, well separated in rapidity from the 
scattered proton, consists of the qq pair. This contribution dominates for the diffractive 
mass M 2 ~ Q 2 . In the large diffractive mass limit, M 2 ^> Q 2 , additional components 
like qqg have to be taken into account. These two components are shown in Figure |6|. 
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The cross section for the diffractive qq production reads M 
da ° L 1 J d 2 rdz\* T>L (r,z)\ 2 a 2 (x,r), 



dt \t=o 16 Ti 

where t is the squared momentum transfer to the diffractive system. Assuming 
exponential dependence on t, e Bot , we divide (|T8| ) by the diffractive slope Bp (taken 
from the experiment) to obtain the total cross section o~tl- The dipole cross section a 
is given by the saturation model (^) which parameters were determined by the inclusive 
structure function analysis. Thus, the description of diffractive DIS (for M 2 ~ Q 2 ) is 
parameter free and differs from the inclusive one by the squared dipole cross section. 

We can perform a qualitative analysis using the formulae from Section [| In the 
DIS case, when 1/Q -C Rq(x), we find for the diffractive qq production from transverse 
photons 



a 



j j °"o °"o °l 




T Q A Rl "~ Q 2 R 2 

r<2/Q 2/Q<r<2R 

By the comparison with (Q), we see that the leading scaling result comes from large 
transverse sizes. The contribution from r < 2/Q is suppressed by the additional power 
of Q 2 (higher twist contribution). Notice that the DGLAP modification of the dipole 
cross section for small r does not influence the leading result. Therefore, DIS diffraction 
is strongly sensitive to the region of the transition to saturation in the dipole cross 
section (^). In fact, the saturated form of a is necessary for the observed in the data 
constant ratio a D /a tot . Considering the leading contribution to (§) and ( |l9|) , we find 

(20, 



o tot \n{Q 2 R 2 {x)Y 

which is a slowly varying function of x and Q 2 . 

Summarizing the qualitative analysis, by exposing the semi-hard (2/Q < r < 2R ) 
and non-perturbative (r > 2R ) regions, DIS diffraction is an ideal process to test parton 
saturation realized in terms of the dipole cross section (|4]). The constant ratio a D /o~ tot 
as a function of x and Q 2 finds a natural explanation in the saturation model. 

In we perform an extensive comparison of the saturation model predictions with 
the data by modelling the diffractive state as shown in Figure |6|. Good agreement is 
found for both the HI and ZEUS data. In particular, in Figure [1] we show the energy 
(xjp) dependence of the diffractive structure function which is determined by the energy 
dependence found in the inclusive structure function analysis. In particular, away of the 
region M 2 -C Q 2 (/3 — > 1) where the higher twist longitudinal qq component dominates, 
we find for the leading twist part of the diffractive structure function 

F 2 D ~ x 1 ^, (21) 

where the "effective pomeron intercept" ctp = l + A/2 « 1.15 (A is defined in Section ||). 
This value is in remarkable agreement with the experimental values: 1.17 by HI [[b 
and 1.13 by ZEUS |5[. 
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We have to add, however, that the treatment of the qqg component, interacting 
with the proton with the same dipole cross section (up to a colour factor) as the qq 
system, goes beyond the saturation model since this component is not present in the 
inclusive analysis. More detailed studies are necessary 

The presented model of DIS diffraction can be extended to include more complicated 
diffractive states, e.g. with partons strongly ordered in transverse momenta (DGLAP 
configuration). In this case initial diffractive parton distributions are directly computed 
from the saturation model and subsequently evolved using the DGLAP evolution 
equations. The results are presented in |I6| . 

8. Conclusions 

The idea of parton saturation realized through the proposed model of the dipole cross 
section turns out to be very successful in the unified description of inclusive and 
diffractive DIS at small x. In particular, the transition to small Q 2 in inclusive DIS, 
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the constant ratio a D /a tot and the energy dependence of diffract ive DIS are naturally 
obtained in the presented approach. From a formal point of view, these ideas allow to 
obey unitarity of the description. The detailed QCD picture of the discussed processes 
is still a matter of intensive theoretical studies. From a phenomenological point of view 



the most promising is the analysis done in [|K], [L7| in which multiple pomeron exchanges 
are responsible for saturation as postulated in our model. In this case the dipole cross 
section (or better the qq - p forward scattering amplitude) obeys a non-linear evolution 
equation. Therefore, the future offers an exciting time for new developments. 
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